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We describe the first DNS-based measurement of the complete mean response of a turbulent 
channel flow to small external disturbances. Space-time impulsive perturbations are applied at one 
channel wall, and the linear response describes their mean effect on the flow field as a function of 
spatial and temporal separations. The turbulent response is shown to differ from the response a 
laminar flow with the turbulent mean velocity profile as base flow. 



Even though a turbulent flow is a nonlinear phe- 
nomenon, a linear response (either in the frequency or 
the time domain) can be defined for it if perturbations 
are small enough. We discuss in this paper the linear re- 
sponse of the velocity vector in the whole flow field and 
in time to small external perturbations applied at one 
wall. The vectorial nature of wall perturbations makes 
the response a tensorial quantity, that we denote with 
Hij to indicate the i-th component of the response to 
an impulsive wall forcing in the j-th direction. Wall- 
based forcing is particularly meaningful if the response 
has to be used in the context of turbulent flow control, 
which is the background of the present study: applying 
the required control at the wall is indeed the easiest con- 
figuration in a realistic flow-control setup. In view of the 
flow-control possibilities offered by modern MEMS tech- 
nology, the linear response function can help considerably 
in controller design, by answering such a basic question 
as which effects are to be felt here and now if a wall 
actuator has been moved there and a given time ago. 

The important role of linear processes in the self- 
sustaining (non-linear) turbulence cycle which takes 
place near the wall has been recently emphasized 1 . Lin- 
ear control has already proven itself succesful in turbu- 
lent flows^. Bewley and coworkers are among the most 
active groups in the field of linear optimal controller de- 
sign for turbulent flows^^. They employ Kalman filters 
and matrix Riccati equations, which require the state 
equations of the system to be known. The (mean) state 
equations are not available, however, for a turbulent flow, 
and must be replaced by a linearized laminar model; in 
addition, the "system noise" is treated as white noise by 
their control design method, and the available informa- 
tion about the actual turbulence statistics does not enter 
the model. Hence these results, encouraging as they are, 
build on an essentially linearized laminar analysis. Re- 
placing the parabolic Poiseuille velocity profile with the 
turbulent mean profile, as done by Hogberg, Bewley and 
Henningson*, certainly improves the effectiveness of the 
controller, but the (mean) effects of turbulent mixing, 
which we aim at describing via the response function, 



still remained unaccessible to controller designers. The 
laminar linearized response has been recently illustrated 
with full detail 6 . The differences between this laminar re- 
sponse and the mean input /output response of an actual 
turbulent flow will be discussed further down here. 

Unfortunately, the obvious definition of an instanta- 
neous linear response is not as useful as might be hoped 
in the context of turbulence: such response is bound to 
exhibit temporal divergence, owing to the chaotic nature 
of the flow. A mean response can however be given a 
precise meaning and measured, either experimentally or 
numerically. One paper which pioneered the approach is 
that by Hussain & Reynolds^, where the response was 
experimentally measured in the frequency domain at a 
given separation and for a few frequencies. 

A few years ago, Quadrio & Luchini 5 proposed a 
method to compute the linear impulse response function 
of a wall-bounded fully-developed turbulent channel flow 
to perturbations applied at one wall. Following that suc- 
cesful proof-of-principle, we proceed here to describe and 
characterize the complete response function tensor, mea- 
sured through a set of purposefully carried out Direct Nu- 
merical Simulations (DNS) of a turbulent channel flow. 
The impact that the availability of such response function 
will have in the field of turbulence control is the subject 
of ongoing work, and it has been preliminarly addressed 
by Luchini, Quadrio & Bewley^, who have been able to 
demonstrate a controller based on Wiener filtering and 
the present response function. 

Let us consider an indefinite plane channel, bounded 
by two walls parallel to the homogeneous directions Xi 
(streamwise) and x 3 (spanwise) and located at X2 = 
and X2 = 2h. The velocity components are u±, 112 and 
U3. To define an impulse response tensor, we input to 
the system an infinitesimal wall velocity perturbation 
Wj (xi , x 3 , t) = Uj (xi , 0, x 3 , t) given by: 

Wj(xi,x 3 ,t) = ej S (x 1 )S(x 3 ) S(t), j = 1,2,3 (1) 

with S denoting Dirac's delta function. The output to be 
measured is the mean effect of this perturbation on the 
velocity field throughout the channel at all subsequent 
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times. The impulse-response tensor Hij so obtained re- 
lates the mean linear response of the turbulent flow to a 
generic input Wj via the convolution: 

Ui(xi,X 2 ,X3,t) = 

J Hij(x\ — x[, X2,x^—x' 3 , t — t')wj(xx, £3, t)dxidx' 3 dt' 

(2) 

Since turbulence fluctuations are large compared to the 
amplitude ej of the external perturbation, which must be 
small enough for the response to be linear and the rela- 
tion |(5J) to be valid, the definition of H cannot be of 
direct use for its actual measurement. However, an en- 
semble average can be used conceptually to define the 
mean response over repeated applications of the impul- 
sive forcing. 

Our first attempt to calculating H has been similar to 
that employed by Hussain & Reynolds 7 , i.e. working in 
the frequency domain: a DNS of a turbulent channel flow 
is performed, where the boundary condition is: 

Wj(xi, X3, t) — ej$m(kiXi)siri(k3X3)sm(u)t), j — 1,2,3. 

Once a suitably small 12 amplitude ej is chosen (which 
in general depends on the forced component), and fre- 
quency lj and wavenumbers fci and k 3 are given, a phase- 
locked average allows the deterministic effect of the per- 
turbation to be separated from the turbulent noise with 
reasonable values of the signal-to-noise ratio S/N. How- 
ever this numerical experiment only yields the response 
function in a single point of the 3d space (ki,k3,u>), 
and we soon realized that the repetition of the computa- 
tion for a number of frequencies and wavenumbers large 
enough to yield a reasonably complete characterization 
of would have been impractical. 

Then we turned our attention to the direct use of (QJ as 
boundary condition, with a suitably small amplitude 
from a computational viewpoint, the ensemble average 
can be replaced by an average over periodic repetitions 
well separated in time, and the complete response func- 
tion is obtained at once. We again realized early that this 
too was going to be an unaffordable simulation: whereas 
impulsive forcing provides in one shot the same amount of 
information as many sinusoidal simulations, it does so at 
the expense of larger nonlinear effects. The correspond- 
ingly smaller allowed ej implies a smaller S /N ratio, and 
the averaging time required to bring S/N within reason- 
able limits becomes unaffordably long. 

We eventually realized that the best of both worlds 
could be obtained by resorting to statistical correlation 
as a method for measuring the impulsive response. It is 
well known from signal theory that, when a white noise 
(i.e. a delta-correlated signal) is passed through a linear 
system, the correlation between input and output is pro- 
portional to the impulse response of the system. We thus 
adopted an indipendently generated random signal as our 
wall forcing, and obtained at once the whole space-time 
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FIG. 1: (Color online) Effect of the white-noise power spectral 
densities ej on the mean friction coefficient Cj, for different 
components of the wall random forcing Wj. C/.o is the fric- 
tion coefficient of the reference flow with no-slip boundary 
conditions. 

dependence of the impulse response by computing such 
a correlation. 

In our method a DNS is performed with a zero mean 
white-noise signal (the output of a random-number gen- 
erator) as boundary condition and the space-time corre- 
lation between this boundary condition (input) and the 
whole flow field (output) is accumulated. Since the ap- 
plied random signal is uncorrelated to the turbulent fluc- 
tuations, the latter will be averaged out just as in phase- 
locking, and the deterministic response will progressively 
emerge while the simulation runs. Moreover, the forc- 
ing power is uniformly distributed over time and space, 
as opposed to what occurs in impulsive forcing, and the 
amplitude can be as large as with sinusoidal forcing. 

The numerical simulations are carried out with our 
DNS pseudo-spectral solver, whose characteristics have 
been described elsewhere^. Of particular relevance here 
is the ability of the code to run in parallel with high ef- 
ficiency. The Reynolds number is Re T — 180 based on 
the friction velocity and half the channel width. The 
domain size is L x = Airh and L z = 4.2h, so that 192 
and 128 Fourier modes (before dealiasing) in the stream- 
wise and spanwise directions respectively, as well as 128 
point in the wall-normal direction, are required to match 
the commonly employed spatial resolution^!.. Peculiar to 
the present simulations is the extremely long integration 
time, about 10 5 viscous time units, required to extract 
the deterministic response from the turbulent noise. This 
time interval is more than one order-of-magnitude larger 
than what is typically employed to obtain converged low- 
order statistics of the turbulent flow. The correlation is 
computed from products in spectral space whenever pos- 
sible; its full y behavior, as well as 81 time separations 
from t + = to t + = 64 are recorded. To minimize disk 
space requirements, a slightly reduced set of 64 stream- 
wise and 84 spanwise Fourier modes is analysed. 

A key step towards the measurement of is the 
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FIG. 2: (Color online) Variation with Xz/h of the absolute- 
value maximum of H\2/t2 in wall-parallel planes. The figure 
refers to i/12 at t + = 30. 




FIG. 3: (Color online) Streamwise derivative d\ip at K2//1 = 
0.1 (shaded contours) of the kinetic potential, given by 
Eqn.([3]l, compared to H12 at t = (line contours). Levels 
at ± 0.02%, 0.01%, 0.005% and 0.0025% of the maximum. 



choice of the amplitude ej of the white noise applied at 
the wall, which must be empirically determined based 
on the requirement that it may yield a linear response. 
A preliminary estimate can be obtained by observing 
whether or not the forcing alters the time-mean value 
of the wall friction. This is shown in figHJ starting from 
the largest ej for which (separately for each forcing com- 
ponent) numerical stability is preserved, we progressively 
halve the amplitude, and observe how the modifications 
of the mean friction induced by the non-homogeneous 
boundary condition becomes negligible: this happens 
only for the smallest values of 6j included in the figure. 

An actual linearity check is shown in figdl where for 
a time delay of t + = 30 the maximum absolute value of 
-^12/^2 m wall-parallel planes is plotted as a function of 
X2 1 h. Linearity requires the curves at different £2 to over- 
lap. This is the case for the two smallest values of e 2 , at 
least in the region of maximum response. At larger x 2 /h 
the curves do not collapse anymore but this is expected, 
since the background noise overwhelms the deterministic 
part of the response and the noise level is different for the 
various simulations (run for the same averaging time) . 

Having assessed linearity, and before turning to illus- 
trate the spatio-temporal behavior of Hij, we devote a 
last preliminary consideration to the response at t = 0. 
The impulse response to wi includes a potential compo- 
nent, that can be computed analytically by solving the 
Laplace equation for the kinetic potential (p(xi, x 2 , X3) 
between two indefinite plane walls. The wall- normal 
derivative 82^ has boundary condition 92^(^1,0,^3) = 
8{x\)5{x3). After Fourier-transforming, the problem sep- 
arates into one-dimensional problems for each pair of 
wavenumbers k\ and k%. Its analytical solution reads: 



Cos1i(k(2 - x 2 )) 
KSinh(2K) 



(3) 



k( 



k 3 . 



where k 2 

In fig[3]the analytical solution ([3]), represented in phys- 
ical space, is compared with H12 measured in the tur- 
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FIG. 4: (Color online) Decay rate of the absolute maximum 
value of H12 : comparison between the full turbulent response, 
the laminar response based on the turbulent mean profile, and 
the laminar response based on the Poiseuille parabolic profile. 



bulcnt flow via the correlation method and shown at 
t = 0: a substantially good quantitative agreement can 
be remarked, except for the lowest contour levels, where 
the residual noise becomes apparent in the turbulent re- 
sponse. Despite the singular nature of the potential com- 
ponent of the response, which manifests itself in a spike in 
the numerically measured correlation, this singular com- 
ponent is faithfully reproduced by our DNS, in which a 
delta-correlated boundary condition is used for w 2 , and 
the correlation between this boundary condition and the 
whole velocity field is accumulated over time. 

We now move on to consider the spatio-temporal struc- 
ture of the full tensor . We shall try to emphasize the 
differences between the computed turbulent response and 
the 'laminar' response employed until now in feedback 
flow-control optimization. To be precise, two kinds of 
linearized response have been previously considered by 
other authors: the solution of linearized Navier-Stokes 
equations about Poiseuille flow and the same solution ob- 
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FIG. 5: Side-by-side comparison of H12 for the laminar re- 
sponse (top) , the pseudo-turbulent response (middle) and the 
full turbulent reponse (bottom). Isosurfaces at ± 0.7 % of the 
maximum, negative values in light gray (enhanced online). 

tained by using the actual turbulent mean velocity profile 
as base flow. This latter pseudo-turbulent response ac- 
counts for the mean turbulent profile but not for turbu- 
lent mixing. While the two responses with the turbulent 
mean profile should be identical at t — and similar at 
very short time delays, they can be expected to progres- 
sively diverge later owing to this difference. That this 
is indeed the case can be appraised from figlH which re- 
ports the temporal decay of the maximum absolute value 
of -H12, the most frequently used component of fly, in 
the whole volume: the true turbulent response clearly 
presents a faster decay rate. From the same figure it 
can be appreciated how the laminar response is differ- 
ent from its companions even at t = 0, and remains 
markedly different from the pseudo-laminar one for the 
considered time interval, thus explaining the control per- 
formance improvement observed in above cited papers 
upon switching from laminar to pseudo-laminar response. 

By observing the three responses in 3d (see fig©, other 
differences can be noticed. In all cases fli2 presents an 
elongated region of near- wall negative u\, but quantita- 
tive differences are considerable. At t + = 30, this region 
turns out to be much longer and narrower for the lami- 



nar response. More importantly, both laminar responses 
- and in particular the one with the mean turbulent pro- 
file - present side regions of positive u\ , that are absent 
in the turbulent one, except for very short time delays. 
From the differences between the true turbulent and the 
pseudo-turbulent responses one is thus led to conclude 
that a possibility exists for further improvements in con- 
trol effectiveness, if the additional information embodied 
in the true linear response can be exploited. 

A comparative glance at all 9 components of fly (fig© 
reveals that they are of significantly different magnitude. 
The components fla turn out to be relatively small, ap- 
proximately 10 times smaller than the others. This par- 
allels a similar observation^ made for the laminar case. 
H2j, i.e. the wall-normal component of the response to 
any forcing, decays much faster than the other compo- 
nents. From a qualitative viewpoint, it thus appears that 
the largest effect with wall-based forcing can be obtained 
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FIG. 6: The nine components of Hij. Top: H\j\ middle: 
i?2j; bottom: Hzj- Isosurfaces of Hn (left) at ± 0.02 % of 
the maximum; isosurfaces of Ha (middle) and His (right) 
at ± 0.2% of the maximum. Negative values in light gray 
(enhanced online). 

with U2 or 1*3 actuation at the wall, and mostly u\ and 
U3 perturbations are introduced into the flow. This qual- 
itative statement is made quantitative by the knowledge 
of the response tensor. 
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Being limited to small control amplitudes is a reasonable reduction 
hypothesis, as long as control is aimed at turbulent drag 



